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Abstract. If the amplitude of primordial gravitational waves is measured in the near-future, what 
could it tell us about bigravity? To address this question, we study massive bigravity theories by 
focusing on a region in parameter space which is safe from known instabilities. Similarly to inves¬ 
tigations on late time constraints, we implicitly assume there is a successful implementation of the 
Vainshtein mechanism which guarantees that standard cosmological evolution is largely unaffected. 
We find that viable bigravity models are subject to far less stringent constraints than massive gravity, 
where there is only one set of (massive) tensor modes. In principle sensitive to the effective graviton 
mass at the time of recombination, we find that in our setup the primordial tensor spectrum is more 
responsive to the dynamics of the massless tensor sector rather than its massive counterpart. We fur¬ 
ther show there are intriguing windows in the parameter space of the theory which could potentially 
induce distinct signatures in the R-modes spectrum. 
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1 Introduction 

It is well known that the pioneering work of Fierz and Pauli [1] sparked a theoretical programme 
aimed at deriving a consistent fully non-linear theory of a massive spin-2 field. This search has 
been refueled by the observation of the current accelerated expansion of the universe [2-4], an 
acceleration which is sometimes predicted in massive gravity models. There were, however, serious 
obstructions to constructing such theories, both theoretical and observational. Of these, the presence 
of the Boulware-Deser ghost [5] and the vDVZ discontinuity which prevented the massless limit from 
reproducing the results of General Relativity (GR) in the regime where they have been verified, are 
the most notorious. 

The ghost-free extension was put forward only very recently [6-8] up to fully non-linear level [9, 
10], and it relies on carefully chosen interactions which are ghost-free by construction and implement 
the so-called Vainshtein screening [11]. Once active, this mechanism guarantees not only the recovery 
of the GR limit, but also improves the stability of the theory under quantum corrections (see, for 
example. Refs. [12-16]). For reviews on massive gravity see Refs. [17, 18]. 

There has been extensive work exploring cosmological solutions in massive gravity [19-35] (and 
its bigravity generalization). If the graviton mass is responsible for the late time acceleration, it 
appears sensible to set it to order of the Hubble parameter today or even smaller. While an ob- 
servationally motivated choice, there are arguments [14, 15] as to why a small mass is a technically 
natural one. Most stringent constraints on the graviton’s mass arise from the physics of solar system 
scales, and therefore late time observations. 

In the case of massive gravity, one may wonder whether there exist further bounds on the 
graviton mass from early time cosmology. On the other hand, in bigravity, which one can think of 
as massive gravity equipped with an additional Einstein-Hilbert piece for the reference metric, the 
dynamics is richer and strongly dependent on the role played by massless as well as massive tensor 
modes. Which observables could be more sensitive to the massive (bi)gravity dynamics? Massive 
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(bi)gravity comes with five(seven) degrees of freedom(d.o.f.) and, in the current data-driven era, one 
might want to test the presence of these additional d.o.f.s at lowest order in perturbation theory. 
However, the dynamics of the scalar and vector sectors can and sometimes is efficiently Vainshtein- 
screened, rendering the cosmological evolution almost unaltered for most part of cosmic history. 
This is precisely what one would want, but ought to verify, if any massive (bi)gravity theory is to be 
employed to describe late time acceleration. 

On the other hand, measuring the amplitude of primordial Gravitational Waves (GWs) can be 
used to constrain various massive (bi)gravity models. This amplitude is traditionally parametrised by 
the tensor-to-scalar ratio, r, which is defined as the ratio of amplitude of tensor to scalar fluctuations. 
Naively, a massive mode decays after crossing the horizon during inflation, lowering the tensor power 
compared to the usual GR prediction. However, there are interesting subtleties related to identifying 
such a signal in the data [36] . In bigravity in particular, the presence of the additional tensor modes 
is not screened in the Vainshtein sense, calling for an investigation of this additional sector. 

In massive gravity there are naturally two metrics (for a general approach see [37]). If only 
one of the metrics is dynamical then there is a clear cut massive eigenstate for the tensors and 
r is sensitive to the corresponding graviton effective mass. In bigravity, however, both metrics 
are dynamical [38, 39], resulting in two sets of coupled tensor modes with time-dependent mass 
eigenstates. These will contribute to the tensor power spectrum in a non-trivial way. It follows 
that the imprint of bigravity on primordial GWs is necessarily less transparent, which has motivated 
further work in bigravity phenomenology in the early universe [40-42] and other contexts [43-53]. 

We focus on a regime of the bigravity theory where several stability requirements are met—see 
Refs. [22, 44, 52, 54], where the background cosmology was extensively studied. In agreement with 
the perturbation analysis discussed in Refs. [55, 56], in that regime the Higuchi bound [49, 52, 57] is 
satisfied and the gradient instability [22] is pushed outside the reach of the effective theory [56]. 

Outline. —This note is organised as follows. In §2 we review massive bigravity in general and specif¬ 
ically the conditions describing the regime of Ref. [56] , focussing on a theoretically viable region in 
parameters space. We choose consistent couplings to matter which fix the way gravity communicates 
with photons that go on to form the Cosmic Microwave Background Radiation (CMBR). In §3 we 
compute the tensor power spectrum and discuss how a measurement of r would impose constraints 
on massive gravity and bigravity. In §4 we elaborate on the possible R-modes imprints of massive 
bigravity in the CMBR signal. We summarize our work in §5. 

2 Review of massive bigravity 

The ghost-free, non-linear extension of the Fierz-Pauli [1] mass term can be written 
TlmGR = f-^b] + ^ n! (4 — n)\ 

where Mg refers to the Planck mass associated with the metric an are constant interaction 
coefficients and the composite tensors JC are dehned by 

JCf^n[9J]=K-X^u With (2.2) 

The scalar interaction potential is routinely and symbolically given in terms of the Levi-Civita 
tensors, Cn[IC] =8 £ KX where the summation of indices is implicit— the interactions are built out 
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of characteristic polynomials of the eigenvalues of 1C, which is at the core of the ghost-free nature of 
the theory. This theory goes generally by the name of massive gravity. 

If the reference metric, f^i,, becomes dynamical, then the theory is referred to as bigravity. 
In this case, the Lagrangian in Eq. (2.1) needs to be augmented by the kinetic term of the metric 
f^lu, given by the corresponding Einstein-Hilbert term (with respective Planck mass Mf), which is, 
to date, the only known ghost-free, Lorentz-invariant kinetic term allowed in four dimensions [58]. 
Therefore, in bigravity models, the Lagrangian becomes^: 

m2 _ 

-^ImBiGrav = ^mGR H- ^[f] ■ (2-3) 

It is important to emphasise the following points about generic theories of massive bigravity. 
First, notice that in these theories there is only one value for the bare graviton mass, m, which 
enters in the interaction potential, C. The two sets of helicity-2 modes generated by the two metrics 
5 ^ 1 / and are associated with a massive and a massless spin 2-fields. The massive set of modes 
has an effective mass, rueg, which is in turn related to m through the an coefficients of the mixing 
interactions Cn- This ruefr is the value which enters e.g. the predictions for observables such as the 
tensor power spectrum and r. Observational constraints are often directly sensitive only to “dressed” 
versions of the graviton mass, a fact we shall make use of repeatedly in this note. 

We also stress that bigravity theories describe two dynamical metrics while being entirely con¬ 
sistent with Weinberg’s theorem [59, 60] which argues that no two massless spin-2 fields mediating 
the long-force gravitational interaction can coexist. In bigravity there is a massive and a massless 
field, so that no violation of such theorem occurs. 

2.1 Bigravity cosmology 

The task of finding viable massive (bi) gravity cosmological solutions has recently received consid¬ 
erable attention, generating a rich and interesting literature [40-42] and other contexts [43-52]. A 
no-go theorem exists for FLRW solutions in massive gravity [34]: there are no FLRW solutions in 
a massive gravity theory with a Minkowski reference metric, Far from being a problem, 

this realization has lead to further work in several different directions. 

Within massive gravity, one may relax the exact homogeneity or isotropy assumption by intro¬ 
ducing some degree of inhomogeneities in the Stiickelberg fields [34], consider open (closed) FLRW 
solutions [21, 25, 61, 62], investigate the cosmology arising from a different choice for / [49, 63] or 
study an extended, ghost-free version of the theory preserving the same five degrees of freedom [64]. 
Naturally, another way out of the no-go is that of considering theories with additional degrees of 
freedom. Massive bigravity and generic multi-metric theories, the quasi-dilaton model [65-69] all 
belong to this class. Alternatively, one may choose a non-trivial coupling to matter [70]. 

In any of the above contexts, once the background solutions have been found, a lot of care 
must be exerted before declaring those solutions viable. Setting aside for a moment observational 
constraints, the study of background and perturbations alone might indeed reveal some pathologies. 
A healthy branch of solutions may support vanishing kinetic terms (the origin of strong-coupling 
issues) or it may lead to perturbations which do not satisfy the unitarity (Higuchi) constraint or 
present a gradient instability. 

Complementing these requirements are those stemming from observational viability: in short 
one would want the “background” cosmic evolution to be essentially ascribable to GR up to the 

^To this one generally adds minimally-coupled matter sector(s). 
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latest era, that of dark energy domination. One should add that, in order to pass, for example, solar 
system tests, an active and effective Vainshtein mechanism also needs to be in place. 

In what follows we will review the analysis and employ the notation of Refs. [55, 56]. Our 
interest lies in the flat, homogeneous and isotropic FLRW universe. We take 


dSg = ff^i/dx^dx^ = —N‘^{t)dt‘^ + a^(t) dx^ 

dsj = fn^dx^dx'^ = —J\f‘^{t)dt‘^ + b^{t) dx^ , 


respectively for each metric. The background Friedmann equations read 


Pm,g + 


m 


Pg 

Mr 

pf 


— PrnJ + j.^2 ’ 

■' K KMf, 


(2.4) 

(2.5) 


( 2 . 6 ) 

(2.7) 


where Hg = a/aN and Hj = b/bM are the Hubble parameters, k = Mj/Mg, and pg and pf are the 
energy density associated with matter coupling directly to the metric 5 ^ 1 / and fgu, respectively. It 
will be convenient for simplicity of notation to define the following dimensionless quantities 

Pm,g — U{^) — - U , PmJ — U (^) , Pm — Pm,g{0 ~Pm,f{C) ) (2-8) 

with primed variables being differentiated with respect to ^ = ^(t) = b(t)/a{t), and 

U{0 = -«o + 4(^ - l)ai - 6(e - lfa 2 + 4(^ - l)3a3 - (^ - 1)V , (2.9) 


where anS correspond to the interaction coefficients, cf. Eq. (2.1). To avoid a number of pathologies 
described in Refs. [22, 56, 71], we choose the so-called healthy branch'^ of cosmological solutions which 
corresponds to the dynamical equation 

Hg = iHf . (2.10) 

It will also be convenient for our subsequent analysis to define 

j{0 = ^ {u{i) - \ V'(o) . (2.11) 

It follows from this definition that J(^) is an implicit function of the an coefficients which set the 
strength of the several interactions in the Lagrangian (2.1). 


2.2 The low-energy regime as a viable region in parameter space 

In bigravity, the mass eigenstates corresponding to the massive modes are generically time depen¬ 
dent [52]. In order to isolate the massless and massive tensor modes at each given time one ought 
to diagonalise the modes using an appropriate basis. It is a convenient feature of the region of the 
parameter space studied in Ref. [56] that the effective mass in this regime is actually constant, which 
makes the analytic derivation of observable quantities, such as r, much more tractable. 

^See e.g. [21, 25, 71] for a more detailed analysis on choosing branches. 
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We shall start by defining the low-energy limit regime (henceforth LEL), investigated in Ref. [56], 
as that satisfying the following condition 


Pg 

m? M| 


<C 1 


and 


epf 

Km? Mg 


< 1. 


( 2 . 12 ) 


Implementing both these relations as well as Eqs. (2.6) and (2.7) is of great consequence for the 
value of the bare mass m and the role of the interactions in Pm,g- The realisation that GR provides 
an accurately verified description for most of cosmic history predating the dark energy era demands 
that Pg dominates the RHS of the Eriedmann equation correspondingly. This, together with the 
regime in Eq. (2.12) translates into the following conditions: 


m Hg , Pm,g < Hg/m ^ Pm,g 1 


(2.13) 


Some of the freedom that comes with having the an coefficients in pm,g is then being put to 
use in taming the “large” value of the bare mass m in Eq. (2.6). We continue below navigating 
the parameter space of the LEL regime to expose its dynamical consistency and, at the same time, 
underline the bounds the latter demands. 

Combining Einstein equations with Eq. (2.10) one derives the following; 


Pm{C) — 


Pg , gV/ 
m‘^M^ Km‘^M^ 


Pm{C) 


< 1 . 

LEL 


(2.14) 


It is intuitively clear then that in the LEL a solution to Eq. (2.14), seen as a dynamical equation for 
is that of a constant ^ satisfying PmiCc) = 0-^ As to the value of ^c, stability considerations 
we will touch upon later on suggest it be order unity. To fully specify the LEL regime, one needs to 
further require the following relation holds true: 


1 > 




rsj 


(Ac) 


Pm,g{Cc) 

Pm,,giA) 


(2.15) 


This can be interpreted as follows: it corresponds to requiring \m'^pmg /<C 1. It holds in view 
of the low-energy Higuchi bound, > 277^), and Eq. (2.6) and can hold in eras preceding dark 

energy domination. Condition (2.15) may be arrived at by enforcing k <C 1 or by appropriately 
choosing the coefficients within pm^g to enforce pm,g ^ p'm,g- Note that the latter option is not 
light on the a^s because of the requirement already in place via Eq. (2.13). 

We note in passing that the lack of gradient instability in the scalar sector corresponds to 
implementing 


c2~l + 


2(c —l)dlnJ 
3 dln^ 


^eiPf+Pg) 

SKM^mlg 


with 


2 _ 2 
rUeff = m 


1 + 


r(0 , 


(2.16) 


and where 




(2.17) 


^If the solution in the low energy regime is not such that = 0, then pm{^c) can always be reabsorbed into a 

cosmological constant contribution. 
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The unitarity bound is automatically satisfied^ upon enforcing Eqs. (2.13) and (2.15). Note the 
appearance of the effective mass, rUes, which is the physically relevant quantity here and corresponds 
to the bare mass dressed by the non-linear interactions which make up the Un potentials. For later 
use we point out here that already at an intuitive level it is clear that a large effective mass would 
give way to a dynamics dominated by the massless modes. One may borrow the familiar -integrating 
out d.o.f.- physical picture to guide this intuition 

As we have seen, working in the LEL regime comes with a simplihed dynamics. It allows for 
example an expansion of all relevant quantities about ^ as a result of which one can write the 
Friedmann equation in Eq. (2.6) as 

= . (2.18) 

and with A = m‘^pm,g{Cc)- We can thus write 

~ _ 

^-M2(i + ^a 

a relation we will put to use in §§3 and 4. In obtaining our estimates for the tensor-to-scalar power 
spectrum, we require the LEL regime to be a good approximation at the very least as early as the 
onset of recombination, when the imprint of bigravity signatures on the CMB occurs. It is easy 
to verify this scale is well within even the most cautious estimate for the strong coupling scale 
^3,f = 

2.3 Coupling massive (bi)gravity to matter 

Before we compute observables, we need to specify the coupling of gravity to matter. When the 
universe was 380,000 years-old, recombination occurred causing the photons to decouple from the 
hot plasma and free-stream up to the present time. The photons of this era are the CMBR we 
observe in the microwave sky and they contain a snapshot of the physics of the early universe. They 
may then encode information about a theory of gravity which is other than GR. 

How does massive gravity couple to the matter sector? Classically, it was shown by Hassan &: 
Rosen [38] that bigravity can couple covariantly to matter without reintroducing a ghostly degree of 
freedom. At the quantum level, one may worry that the irrelevant interactions which make up the 
non-linear theory might not be under control.® 

In bigravity or multigravity theories, one could wonder whether all the dynamical metrics could 
democratically and covariantly couple to the same matter sector. Unfortunately this induces a ghost 
at an unacceptably low scale [70, 76], making the model unstable. As a result, either each dynamical 
metric couples to its own matter sector, or the coupling of gravity to matter is made through a 
composite metric. In the last case, however, FLRW ansatzs for both g^i, and require matter to 
which they couple to be identified as a dark sector, [70, 77] which can later on couple to standard 
matter. This would require knowledge of this additional coupling. 

For simplicity, in what follows we couple matter to one of the metrics, 5 ^,^. 

^See Refs. [49, 52, 57, 72] for more details. 

®See also [83] on this point. 

®In general, the quantum stability of field theories which rely on large derivative self-interactions, as is the case 
of massive gravity, is not trivial. However, recent work has shed some light on the role of the quantum mechanical 
realisation of the Vainshtein mechanism in protecting the theory—see Refs. [12, 13, 16, 73-75]. 


(2.19) 
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The matter Lagrangian is fully described by 


<5matter — j d X \J /Imatter(fl') 


( 2 . 20 ) 


where 4>i labels the matter field species, including photons which encode the physics of the early 
universe after recombination occurs. 


3 Primordial gravitational waves 

Primordial gravitational waves are seeded by tensor modes of the primordial perturbation. In massive 
gravity, there is only one family of modes, which we shall denote by h^y whereas for bigravity there 
are two sets of modes at play, h^y and ifj,y. Therefore, the fluctuations of both metrics, 5 ^ 1 / and /^i/, 
now contribute to the power spectrum of primordial GWs. 

We start by reviewing the cosmological perturbation analysis performed by De Felice et al. [56]. 
We write the metric in ADM variables 


g ^y = — dt^ + a{t)‘^ {'^ij + hij)dx^ dx^ and fij = —N‘^dt^ + + (.ij)dx^dx^ (3.1) 

where ^ij is spatially flat, a and b are scale factors and N and M are lapse functions corresponding 
to each space-time metric. The quadratic action of the tensor modes becomes 


5, 


( 2 ) 


tensors 


Mo f A 

/ d^xiVa^^ 


iV2 




- hij + 




1-3 


AA2 62 


+ t 




-m2r(0 


(3.2) 


where 7 can be set to unity, denotes the spatial Laplacian and 

Ma 

Dotted quantities are differentiated with respect to cosmic time. The action above (3.2) features 
non-trivial interactions between h and i fluctuations, which makes analytical predictions far from 
straightforward. To make progress, it is be useful to find a regime in which a diagonalisation of this 
action is possible, as argued before. 

Considering only the transverse and traceless perturbations as the relevant degrees of freedom, 
it is convenient to choose a different basis for the tensor perturbations as follows [56]: 

Hr. = hij - lij and H± = ■ (3.4) 

We stress that the validity of this diagonalization is limited to the LET regime. The mass eigen¬ 
states of a bigravity theory are in general time dependent and so will be the diagonalization basis. 
Nevertheless, in what follows we will employ the customary inflationary normalization for the tensor 
wave functions whose behaviour at recombination is what determines possible CMBR imprints. 

Before proceeding with the analysis of tensor modes in the LET configuration we pause here to 
stress that any instability in the tensor sector should not go unnoticed and it is different in nature 


(3.3) 
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from the one concerning the other d.o.f.s. Indeed, no matter how efficient the screening of scalars 
and vectors, the two copies of the tensors at hand can be safely considered unscreened. The act of 
normalizing things “as usual” for inflation and considering the tensors diagonalized from the onset 
of the LEL may be insensitive to part of the dynamics which goes beyond the LET regime', as is 
exemplihed by the work in [40, 78], which respectively focused on different branches of solutions. 

Working at lowest order in perturbations, the simplified action for the tensor modes becomes 


c.(2) 

*^tensors 


- 




Mi \ HiHt, + ^V^R, 


V 


+ m ! U'lHr + ^ 


(3.5) 


As we shall see shortly, it is this mass parameter nies associated with the tensor modes which 
enters observable quantities and which can be used to constrain the parameters of the theory. The 
normalisation of the kinetic terms is set by the mass scales M+ and M_ given by 

Ml = {l + Ke,)M^ and Mi = . (3.6) 

In perturbation theory the matter Lagrangian (2.20) is given, to lowest order, by 

<55^atter = j . (3.7) 


The usefulness of the {+, —} basis is clear: is the massless mode whereas the massive mode 

is R~-. In light of our choice of coupling to matter (3.7), photons do not couple to a massive nor 
massless graviton. Instead, they couple to a linear combination of massive and massless modes in the 
diagonal basis described in Eq. (3.5). Importantly, the massive mode is associated with a redressed 
mass meff, which albeit related, is not the graviton mass. 

In terms of the diagonalised variables of Eq. (3.4), the coupling to the matter sector is 

= I X' . (3.8) 

This means that gravity and the history of the universe will, in principle, change. To guarantee that 
the cosmology remains invariant so as to agree with ACDM and to reproduce the Newtonian limit, 
we need to require that the Planck mass (or equivalently the Newton’s gravitational constant) is 
essentially the same as measured in solar system scales. Consider two mass tests described by the 
energy-momentum tensor above and subject to gravity set by this matter coupling. One may easily 
derive that the effective Planck mass associated with this modified gravitational force is actually the 
Planck mass associated with the metric 

AIpi,efr = Mg . (3.9) 

^This is important because the coupling among the tensor modes, although very small, in time can have significant 
effect once higher order in perturbations are taken into account. This again points to the importance of the initial 
value problem as emphasized in [40]. See also [42] for a step in this direction. Note that the scales probed by [42] are 
far beyond the reach of the LEL regime. 







So far, all we have assumed within the LEL was that ~ 0(1)) but there were no stringent 
constraints on the value allowed for n. If we demand the cosmological evolution to be the same for 
most of cosmic history, as well as negligible modifications to the Newton’s gravitational constant, 
then from the Friedmann equation (2.19), we ought to require 

K <C ~ 0(1) . (3.10) 

Considering the uncertainty associated with the empirical determination of the value of the 
Newton’s gravitational constant, we estimate k to be 1 part in 100,000. One might be worried by 
the realization that this regime corresponds to a small Mj <C Mpi. This results, in turn, in the 
lowering of the naive strong coupling scale = Asj.® This worry is legitimate and we just 

note here that the scale we are negotiating with is that of Hr- One can easily check that even if the 
strong coupling scale is as low as A 3 J, the theory is still predictive at the time of recombination, for 
which the graviton mass can be chosen such that Hj- <C A 3 J, a condition which is automatic in the 
LEL regime. 

Other d.o.f. — Before proceeding with the study of observables related to tensor perturbations we 
pause here to comment on the role played by the other degrees of freedom in our setup. Clearly 
massive gravity alone already spans 5 d.o.f. whose effect on cosmological evolution up to recombina¬ 
tion must be accounted for. So far we have largely neglected the dynamics of the scalar (helicity-0) 
and vector (helicity-1) sector. Their dynamics has already been the subject of careful investigations 
at the lowest orders in perturbation theory hinting at early(late) time instabilities. However, such 
an analysis may or may not capture the full physical picture, especially if limited to low orders in 
perturbation theory. It has further been suggested that a way around these issues may be found 
in a restricted pool of favourable initial conditions and that the so-called initial value problem for 
bigravity needs further study [40] (see also Ref. [81] for related work). We can only but agree on 
this latter point. A step in this direction was taken in [42] where initial conditions are discussed in 
an inflationary context. This implicitly assumes that the effective strong coupling scale is as high as 
Hint- 

In this manuscript we take the view that, as is often the case in solar system dynamics, an 
efficient® Vainshtein mechanism will make use of the non-linearities of the theory to milden the role 
played by the scalar (vector) sector in the cosmological evolution, effectively screening them. Although 
suggestive results are present in the literature, at this level this is indeed just an assumption we make 
here, and not a small one at that. 

We now turn to the observable spectrum arising in this bigravity model. The coupling in 
Eq. (3.8) dictates the power transmitted by the tensor modes (both from the massive and the 
massless spin-2 fields) to the CMBR photons. Before presenting the formula for r, we make a short 
digression into the power spectra of massive and massless tensor modes which would individually be 
imprinted in the CMBR. For bigravity theories, the power spectrum will be a composite measure of 
both signals. For the purposes of our estimates, it suffices to work at lowest order in the slow-roll 
approximation, which is assumed to hold throughout inflation. 

Power spectrum of massless modes. —The analysis of the perturbation theory for massless tensor 
modes is well known in the literature [82]. In bigravity theories, the massless mode is represented 

®Here we are focusing on the case when the strong coupling scale can be as low as {m^. Being a very low 
scale, this could prove difficult for the phenomenology at and above that energy scale. 

® Incident ally, we note that an active screening in the LEL demands [55] that |dlnJ(C/dln(C| 1, see Table 1. 
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by in Eq. (3.5) and their evolution is fixed by the time they cross the horizon, after which their 
amplitude remains constant thereafter. It is therefore sufficient to determine their super-horizon 
evolution. On super-horizon scales, the resulting power spectrum is given by 

^massless = i/+(k 2 )) = (3.11) 

Here we adopt the usual notation that starred quantities are evaluated at the time of horizon crossing 
for the mode k, and we use H = Hg to avoid clutter. Notice the different normalisation arising from 
Eq.(3.5) of the tensor spectrum from the GR one, which we have denoted by E’t,gR) where the only 
scale is H. 

Power spectrum of massive modes. —Massive tensor modes have their evolution determined by 
their corresponding mass horizon, and they do not get frozen on super-horizon scales but rather 
oscillate. These modes are denoted by H~ in Eq. (3.5) and give the following power spectrum 

Rmassive = (^jj(kl) H-i (k2)) = ^ -ft,mGR (3.12) 

where Pt,mGR stands for the power spectrum obtained from the standard GR action for tensors 
where there is an additional mass term (in this case rueff). Notice again the different non-trivial 
normalisation of the power spectrum (arising from Eq.(3.5)). 

3.1 Tensor-to-scalar ratio in bigravity 

The power transmitted by the tensor modes in this class of theories is governed by how they couple 
to matter via Eq. (3.8). It follows that the primordial power spectrum of tensors is 

(h,,(k0^,(k2)) = (F+(ki)if+(k 2 )) + (^i^(kl)F^(k 2 )) (3.13) 

with the individual power spectra given by Eqs. (3.11) and (3.12). The fact that this prediction is 
different from that of GR is only one source of modification to the tensor-to-scalar ratio, which is 
also sensitive to the scalar perturbations. Since the Eriedmann equation is still slightly changed, this 
will reflect on the matter power spectrum and therefore on the primordial scalar power spectrum. 
In particular, from Eq. (2.19) the energy density in matter is 

^matter = ^ ^ ^2 ’ (3-14) 

which results in the scalar power spectrum changing to 

^’s.biGrav = . (3-15) 

l-hAt^c 

where the subscript s refers to scalar. Gonsequently, these theories predict 

^'BiGrav = J'GR + («^?c) ^mGR j (3.16) 

where rcR denotes the tensor-to-scalar ratio in GR (that is, assuming a massless graviton) whereas 
^’mGR corresponds to that in a theory of a massive graviton. Under the assumption that the cosmology 
is not signihcantly changed and in the limit where <C 1 , it follows that the prediction for r in 
these theories is the same as the original for GR. 


- 10 







3.2 Summary on bounds in parameter space 

To obtain the observable predictions in bigravity theories, our starting point was the region in 
parameter space in the theory where known instabilities are absent [56] which we have further 
restricted due to the k <C 1 requirement. Our results reflect this choice. We summarise on Table 1 
the parameter space available so that the theory is free of known instabilities. This entails that, even 
though there is in principle more freedom to choose from the higher number of degrees of freedom 
in bigravity then in massive gravity, absence of instabilities plays a rather restrictive role. 


Stability requirement 

Parameters 

Bound 

Higuchi bound 

meS 

mis > 

Absence of strongly coupled perturbations 

Hg 


No gradient instabilities & active Vainshtein 

m, H 

m'> H, dlnJ(^)/dln(^) » 1 

Standard Friedmann evolution 

Mg 

K 1 with ^ ~ 0(1) 


Table 1. Collection of stability bounds for the phenomenology discussed in this paper. Notice that these 
conditions correspond to relatively strict bounds on the a„ interaction coefficients in the massive bigravity 
Lagrangian. This can be easily seen by expressing J in terms of via Eqs. (2.9), (2.11). 

The tensor-to-scalar ratio is not the only observable which can provide constraints to this class 
of theories. Interesting imprints in the CMBR are also forecast which, provided they affect observable 
modes, can reveal interesting physics about the early universe. We turn next to such characteristic 
signatures in the R-modes signal, which were originally unveiled in Ref. [36]. We adapt their results 
to discuss their implication on the bigravity models we discussed here. 

4 Possible Signatures in the CMBR spectrum 

The study of the CMB spectrum plays a fundamental, unparalleled role in modern cosmology. It 
essentially represents an open window over the past of the Universe, taking one all the way back 
to the recombination epoch. It is then clear that, in the study of how a massive (massless) graviton 
may affect the CMB radiation, the crucial scale one will have to negotiate with is the size of the 
Hubble radius at recombination, As we shall see below, a simple description emerges for 

the primordial tensor sector contributions to B-modes: one can think of two main regimes for the 
effective tensor mass m-efr, with the transition value for TUeff being itself. 

In deriving the results for this section we will heavily rely on the work in [36]. In there, the 
authors present a phenomenological approach aimed at placing bounds on the graviton mass which 

^°This condition specifies the so-called healthy branch. 

^^This last assumption is to be understood within the LEL [56]. 

^^The contribution from reionisation is relevant for small £. In the massless case ^ > 20 is sufficient to ignore this 
contribution. The numerical analysis of [36] supports the approximation scheme neglecting the reionization contribution 
also for a large range of effective mass values. 
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can be adapted to our case. Their results are quite general and hold true for the tensor sector as 
long as the tensor wavefunction satisfies the usual tensor e.o.m., just as the one that follows from 
Eq. (3.5). As per Section 2, the model under scrutiny here is a bigravity theory, it counts seven 
degrees of freedom which include massive and massless tensor modes. 

As documented in Eq.(3.4), both the massless and the massive tensor modes, decoupled from 
one another in the basis, directly couple with matter. The crucial realization in what follows 
is that there exist important regions in the parameters space where, to a good approximation, the 
result of the two contributions is additive not just at the level of the action^^, but it remains so 
up the source term and therefore eventually propagates all the way to the expression for the 
coefficients ^ . We will argue in particular that this is the case in the mass regime which would 
in principle generate the most prominent effect in the CMB. 

As we have seen, it turns out that in the LET, a regime we chose in order to avoid known 
instabilities, the relative coefficient regulating the contribution of the massive tensor modes (as 
opposed to the massless ones) is very small because it depends linearly on k <C 1. This will result 
in a very hard-to-detect massive tensor sector. We will see below that in a specific mass range a 
massive graviton actually enhances the gravitational signal by two orders of magnitude as compared 
to the massless case. However, in our case this enhancement is no match for the oc k suppression. We 
nevertheless provide a detailed discussion also in the hope that it will be later applied to a scenario 
where the relative coefficient is order one so that “massive” imprints would be more conspicuous. 

4.1 B-modes 

A clear qualitative understanding of the dynamics of tensor modes may be arrived at by classifying 
them (see Fig.l and [36]) according to whether or not and to where these modes are relativistic (in 
the sense of satisfying q^/a(r)^ 3> 

• Modes relativistic at recombination belong to Region I 

• Modes which are non-relativistic already early on, before entering the horizon, reside in Region 
II 

• Momenta entering the horizon as relativistic but turning non-relativistic by the recombination 
epoch populate Region III 

This pictorial view refers of course to the massive modes but in the additive regime one may reason¬ 
ably expect that the exactly massless modes will generate the usual imprints to be superimposed, 
weighted by a relative coefficient, to the massive modes signatures. 


Note that the presence of the fractional relative coefficient in front of the massive contribution to 
e.g. in Eq. (3.8), cannot alter the description in Fig.l Indeed, the latter is obtained by 
judiciously comparing among each other the value over time of q/a{t), H{t), mes■ These quantities 
squared all appear in the massive tensor modes equation of motion and their relative strength, 
unaffected by the fractional coefficient, signals the regime one is working at (i.e. (non)relativistic, 
inside/outside the horizon). 

^^This is trivially true as we linearly couple gravity to matter in Eq. (3.7). 

^“^What it does alter is the relative weight of the massive tensor sector in determining the gw signal. 
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0 100 200 300 400 500 

I 


Figure 1. Above is the effective mass Weff in units of the Hubble rate today, Hq, versus multipole momenta 
The upper limit of Region I is obtained by requiring that the multipole for which i ~ ( 7 (to — t^), 
coincides with ^O) that is the multipole corresponding to a physical momentum the size of the mass Weff at 
recombination. The border between Region II and Region III is obtained by identifying the multipole ^rn 
corresponding to the momentum that becomes non-relativistic at the same time it enters the horizon [36]. 
Two horizontal lines show qualitatively different regimes for the effective mass: m> stands for mass larger 
than Hj. and, complementarily, . 


Let us now discuss how one may generate Fig.l and then the consequences it can possibly entail 
for bigravity CMB signatures. Region I denotes modes relativistic at recombination. For a generic 
multipole moment one has, in conformal time, that I ~ g(ro — Tr). The role of the (tq — t^) factor is 
clearly that of accounting for the time evolution to the present day. In our setup this evolution can 
in principle depart from that of ACDM. 

Working in the LEL regime and implementing Eq.(2.I9) (with k <C 1 and Agfj ~ A = m?pm,g[^c] 
so that the form of the Eriedmann equation reduces, formally, to the standard one) then, in multipole 
language the non-relativistic threshold is reached [36] at recombination whenever: 


£ < meffa(rr)(ro 


mes 


da 


3.3- 


rries 


Hq[ 1 Zr) J + fir Ffo(l + ^r) 


(4.1) 


In the numerical calculation we have set flAes ~ ~ 0.73 because, as we have seen, the 

background evolution of the LEL regime mimics that of ACDM in our setup. 

It is also worth pointing out that the result of a (4.1) without a c.c. term would amount to changing 
the numerical factor in front from 3.3 to 3.7, a mere 10%, reflecting the fact that we owe most of the 
evolution since recombination to the content in . 

A recombination redshift of Zr — 1088 delivers the yellow shaded area in Eig.(l). Note that for 
very small masses, say nieg < 10^Hq, the value of £o will be so small that the B-mode spectrum will 
be unaffected by non-relativistic modes Incidentally, for mefr = Hr one finds that £o ~ 64, which 
is the only point shared by all three regions in Eig.l. 

^®This is true exactly only if the contribution at reionization is neglected. The fact that non-relativistic modes play 
no active role is intuitively clear upon noticing that the mass rrieff ~ 300 i4o corresponds to the size of the visible 
universe at recombination. 
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The first, most straightforward, realization is that modes which are relativistic at recombination 
are entirely insensitive to the presence of the effective mass nies- On the other hand, one can always 
probe the non-relativistic regime through modes satisfying q/a{t < t^) < rUes. 

In the m< range there are only two possible configurations: modes relativistic at recombination 
(Region I), or non-relativistic momenta which only re-enter the horizon after recombination (Region 
II). In this mass range modes are indeed forbidden from stepping in the horizon as relativistic and 
slowing down to non relativistic before recombination. Indeed, consider a q/a{t < R) > nies', this 
mode can easily transition into non-relativistic after some time. On the other hand, because H(t) 
has a steeper time dependence than a~^{t) and becomes H = > m,efr at recombination, it must 

be that if the mode q in question is outside the horizon it will stay out past the time it becomes non- 
relativistic. This dynamics then fits the description of modes populating Region II. The signatures of 
massive tensor modes in the m< range become distinct from those of the massless case only for long, 
outside the horizon, wavelengths. The source term which feeds the expression for the multipole 
coefficients depends on the time derivative of the primordial tensor perturbations hij, which is in turn 
related to the term q^ + rngfr^ a(tr)^. Clearly, at very large wavelength (from the onset of what we 
call the non-relativistic regime) it is the mass contribution to dominate and provide an enhancement 
with respect to the massless case. This is precisely the low-I" plateau found in [36]. 

For completeness we report that, as one raises the value of the effective mass, a new qualitatively 
different possibility emerges: relativistic modes might enter the horizon and slow down to become 
non-relativistic before recombination. This is the dynamics which characterizes Region III. We refer 
the reader to Ref. [36] for details on how to derive the border between Regions II and III. Our interest 
is focused on the modes of Region II which, as we have briefly reviewed, contribute to a plateau in 
the CMB tensor spectrum. 

Having seen how Fig.l provides an understanding on the dynamics of massive tensor modes, we 
pause here to note that this understanding is necessarily qualitative in nature and serves its purpose 
quantitatively only in the asymptotics. We have shown in detail below Eq.(4.1) that a very small 
effective mass, below 10^Hq, will witness most modes being relativistic and, as for the non-relativistic 
ones, those will not correspond to a high enough i so as to leave any marks on the B-modes CMB. 
It is safe to say that such a small effective mass would not be detected 

Navigating Region II for larger meff will eventually lead to the “large wavelength outside-the- 
horizon” enhancement mentioned above, an effect propagating all the way to the B-modes spectrum 
and shaping up as a the low-£ plateau in the ^ multipole coefficients (for a fixed nies of this 
size and higher is one would step into Region I). Determining exactly the onset of this enhancement 
is a task best performed through the use of software such as CAMB [? ], see [36]. The result is 
that a plateau starts emerging at about meg ~ 1.2 x its effect being most striking (two 

orders of magnitude larger than the standard massless tensor signal of GR) at mes ~ 1.5 x lO^F^o 
only to weaken and eventually become suppressed with respect to the massless signal as soon as 
rUeS ~ 3 X lO^FIo- 

We stress here that the plateau is a clear-cut effect that, if detected, would represent the most 
prominent CMBR signature of a massive theory of gravity. Most importantly for our analysis, the 
fact that this effect generically takes place in the low-.^ regime will, as we shall see, guarantee that a 
bigravity model may also lead to such an imprint in the same effective mass range. Another thing to 
keep in mind is that in our case, as opposed to the analysis in [36], the unitarity bound sets a strong 

^®This statement is all the more appropriate in a bigravity setup where the massive tensor modes do not have the 
full weight they enjoy in a purely massive gravity theory . 
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upper bound on rrieff, of the order H^. This requirement stems from the helicity-0 mode analysis and 
is therefore not necessarily present in e.g. Lorentz breaking theories of massive gravity that inspired 
the work in [36]. 

Proceeding with the analysis at larger effective mass values one will see a suppression of the 
signal. The reason for the asymptotic suppression at meg S> is that non-relativistic modes will 
start oscillating sooner and sooner outside the horizon with increased frequency rngg leading to an 
averaging-out which amounts to a strong suppression of the signal. Determining exactly by what 
value of meg this effect will take over the enhancement is a task beyond the scope of the present 
work. 

The additive regime and the bigravity signal. — The reasons Region II is of particular interest 
for us are manifold: besides being responsible for an intriguing \ow-^ plateau in massive gravity, this 
is an area in (( 7 , meg)-space whose contribution to multipole coefficients can be treated as “additive” 
to a good approximation in the case of bigravity. 

As mentioned above, where we part ways with the work in [36] is in considering two sets H~^^~ 
of tensor modes. As a consequence, our source term will consist also of cross terms. Crucially, in 
a low-g and outside-the-horizon range such as Region II, the contribution of the massive tensor modes 
far surpasses (assuming the contributions are equally weighted) that of its massless counterpart. As 
a matter of principle then, not just massive gravity, but also bigravity can lead to distinct imprints 
in the B-modes spectrum for the appropriate meg range. 

In bigravity though, the nature of the signal will also depend on the value of the relative 
coefficient between the massive and massless tensor modes coupling to matter via e.g. Eq.(3.4). 
Schematically: 


Signal 


('^source)lm=0 + Crei 


1 2 


source/ I^efF 


(4.2) 


The first thing to keep in mind is that the unitarity bound immediately sets the dynamics in 
the meg > Hr region. A Orel S> 1 would produce a signal profile with a very substantial overlap with 
the massive gravity one, with a low-.^ enhancement at the specified mgg range and a suppression in 
the asymptotics of a very large meg. For Orel ~ 1 and down to values satisfying Orel ^1/10 one 
would still be able to see the low multiple enhancement for the appropriate meg; this is because the 
maximum signal enhancement due to a mass (when mgg ~ 3 x 10 ®/Fq) is two orders of magnitude 
larger than the would-be massless signal. As for the mgg > 20Hr range in this configuration some 
suppression is to be expected but, again, this regime is best understood by running the software^'. 
Finally, the Crei ^ 1/100 configuration is bound to generate a signal with almost complete overlap 
with the GR profile. 

In general then, the difference between a bigravity theory and its massive gravity limit would 
most clearly manifest itself as a lesser enhancement in the low-£ plateau for the former and cor¬ 
respondingly a less dramatic suppression in the large mgg region of the bigravity parameter space. 
From the perspective of CMBR signatures one may summarize these findings as evidence that bigrav¬ 
ity theories generate imprints which are overall less sharp than those originating from pure massive 

anticipate however that in the rrieff S> Hr regime the source function might well receive an important 
contribution from what we call the -cross terms- of the two tensor sectors and therefore one may not rely on the 
massless/massive modes additivity any longer. On the other hand, it is important to stress that a more detailed 
analysis of the dynamics for which cross terms play a leading role (a purely bigravity effect this one) is bound to lead 
to an interesting characterization of further signatures of bigravity theories. We leave this to future work. 
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gravity. This is justified already at an intuitive level because the bigravity theories space spans cor¬ 
ners of pure massive gravity (e.g. the Mj —)■ oo limit) but is also endowed with two additional degrees 
of freedom which account for the dynamics of the additional massless tensor modes. In principle, 
depending on the relative coefficient Crei of the modes contribution one might make bigravity 

imprints as sharp as those for massive gravity; on the other hand, those signatures can be “watered 
down” towards the purely massless spectrum by a different judicious use of the same coefficient. 

Zooming in the LEL regime, one can see that Crei ~ Cl(l) • k < 10“® thus concluding that the 
imprints of bigravity in the LEL are very hard to probe and the CMBR signal is expected to mimic 
GR. What led to our restricted parameter space were, in addition to our working in a low-energy 
regime, the requirement of an almost ACDM background evolution combined with the constraints 
on the effective Planck mass. In turn, this resulted in a very small relative weight for the massive 
tensor modes at recombination. 

Our setup should by no means be thought of as the only cosmologically viable option. Much 
more work is needed in this direction. It is indeed quite possible that another stable region may 
be found in the future whose domain includes a massive tensor sector generating more prominent 
signature in the CMB. 

5 Summary 

The observed accelerated expansion of the universe has reignited the research aimed at finding 
compelling theories modifying GR in the IR. Theories of massive gravity are a sure candidate for the 
part. Their phenomenology is currently under intense scrutiny. Most studies so far have focused on 
the late-time cosmological dynamics in these theories. In this paper we took a different approach and 
asked how early-time cosmology dynamics and constraints would reflect on the bigravity parameter 
space. 

Our starting point has been a theory of bigravity, where both metrics are taken to be ELRW 
and matter only couples to one of them. We studied the predictions for the tensor-to-scalar ratio, 
r, in bigravity in a specific low-energy limit. In this regime, the two copies of tensor modes, which 
are generally coupled, can be diagonalized in a time-independent fashion and solved for. It is then 
possible to write r as a linear combination of contributions from the massless and massive modes. We 
find that in this regime, once additional constraints are imposed, the contribution from the massive 
tensor modes is suppressed. 

We further showed that, although the massive sector of the theory can in general leave quite a 
distinctive imprint on the B-modes profile, the region in the parameter space we have been probing 
supports a suppression of the massive tensor modes contribution to the overall signal in favour of a 
GR-like profile. 

We are led to conclude that a very efficient Vainshtein screening under specific^® conditions 
leads to a hard-to-detect, as far as GMB data is concerned, bigravity imprint. One could say that 
the massive tensor sector would bypass, rather than pass, the cosmic background radiation test by 
means of a small k. A posteriori this is not surprising: the assumed strong Vainshtein screens the 
non-tensorial d.o.f.’s and k <C 1 effectively reduces Hf-, the massless modes, to hij and suppresses 
H~-. This is morally very close to GR. On the other hand, the k <C 1 condition stems only from late 
time constraints on the effective Planck mass and what looks like GR at recombination may reveal 

Although, as pointed out in [40], an in depth study of the initial data problem is needed in order to decide what 
really is tuned and what is not. 
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itself more at solar systems scales where an analysis of the screening mechanism is more accessible 
via the decoupling limit. 

Further investigation is required to ascertain the degree to which an active Vainshtein can 
screen at the decoupling limit and, possibly, away from it. The latter task is especially interesting 
and important for the analysis presented here, but also especially complicated. Indeed, in our setup 
a small Mj is expected and it corresponds to a small naive strong coupling scale A^j <C A 3 . We 
leave this to future work. 

Relation to recent works. —The phenomenology of bigravity theories has seen increased recent 
interest. Some works have some overlap with our own. Comelli et al. [71], Koennig et ah [80], 
Lagos et al [40], Cusin et al. [78] and Amendola et al [41] employed linear perturbation theory to 
understand what pathologies may arise in each sector of bigravity. These works vary depending on 
the branch and region they focus on but overall they span several branches of solutions including 
the one employed here. A common conclusion is that further investigation is needed on the initial 
value problem, as already stressed in [40]. A step in this direction is represented by the work in [42] 
by Johnson et al. In particular, [42] identified a background which generates GR-like results under 
a large pool of initial conditions during inflation. 

The low-energy regime (LEL), employed throughout this manuscript was introduced in de Felice 
et al. [56]. To the parameter space therein we superimposed additional constraints (e.g. k < 10“^) 
that lead to a reduced space where the massive tensor sector is suppressed. In this sense then, we 
naturally make contact with the work by Akrami et al. in [83]. 

Our working in a restricted LEL regime means the scale at which one imagines initial conditions 
are set is far beyond the reach of our approximation and above our setup strong coupling scale so 
that our contact with the results in [42] is an indirect one. 

While this work was nearing completion, a preprint by Sakakihara et al [84] appeared in which 
the tensor power spectrum was derived in a specific bigravity model, corresponding, in the notation 
of Eq. (2.1), to an ^2 = 0. The Authors considered the so-called healthy branch of solutions and, in 
as much as there is overlap with our work, they reached conclusions not unlike ours. 

A number of earlier works also focused on massive gravity signatures in the CMBR. Ref. [36], as 
we have seen, pointed out the existence of a B-modes enhancement in the low multipoles for a specific 
mass range. The study in Ref. [28] concerned a general analysis of the profile and detectability of 
the gravitational wave signal arising from a time-dependent mass term in massive gravity. 
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